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We present a study of the Adomian's Decomposition Method (ADM) applied 
p I ■ to the Hamilton- J acohi equations ut + H {ux) = 0. 

, We recall the well known characteristics methods in the case of this type of 

equations to justify the existence or not of solutions. This yields that the ADM 
gives efficient solutions in time only in ]0, T*[, where T* is the critical time 
of our equation. 
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1 A presentation of the ADM 
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^ ' We consider in this work the Hamilton- Jacobi equation 



\ Ut + H{ux) =0 xeR, t>0 . . 

" \ u{0,x) = uo{x) X ^ ^ 

where u = u{x,t) and Uq = Uo{x) are functions on a Hilbert spaces, and H a non 
linear operator. 

d 

Let take L = — as the linear operator. The principle of the ADM relies on the 
ot 

assumption that the solution u of equation ([1]) can be set as a serie u = Uk. In 

k=0 

[1] the authors considered the non linear part as a function of Ux. This yields to: 

\dx J ^ y dx ' dx J ' 

and 

Aq = H i j with Mo = "^(0, x) 



= E(/c + iK+i^A„ (2) 

fc=o "^fe 



where u], = — — . 

" dx 



Several authors give some formulas on the Adomian's polynomials A^. in the case of 
an operator depending exclusively on u (see for example [21 II])- For more related 
works in this field one can see |3J, [5J, [8], [Tl], [15] and [T6] . 
We present here other formulas more suited to our analysis. 

Theorem 1 Let's consider that the operator H of ([Ijj is indefinitely differentiate 
in Frechet sense . 
We have: 

Ao= H (u'o) andA^ = Y,— J2 ^^^^ («o) > «P2 ' 

k=l pi+P2+.-.+Pfc=n 

dUk 

dx 

Proof: 

" dv- 

We note V'(A) = E ^'tt • Thus 
i=i 

n\ dX"" n! rfA" ^ ^ 

n 2 

but [H o ^](A) = [H o ^](Ao) + E 7t[^ ° (Ao) (A - Ao, A - Aq, , A - Aq) 

k=i - ' 



where H^^^ [u'q] is the k-th derivative at w^, and 



k times 
then we deduce that: 

H{^{X)) = H{^{X,)) 

n 1 1 
k=l^- P1+V2+ ■■■+Pk=n Pl-P2----Pk- 

For Ao = 0, we have: 



^1 A;! +p,=n Pi'P2! •••Pfc! \dxj\ dx dx 



Vk 



H 



dx 



fc! Pi!p2!...Pfc! \dx)\dx'""' dx J 

duo" 



H 



dx 



t'l p,+p,^+p,=n \dx)\dx'""' dx 



Thus: ^Hm)) = n! E A E ^^^^ --^ ^ 

.=1 p,+p./::'+p.=n \dx)\dx' ' dx 



This ends our demonstration. 
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Theorem 2 For the problem P^, we have A^-x = Un{x) 



j.n— 1 



(n-1)! 



and Un(x,t) = Unix) — Vn > 1, with 

n'. 



n-l 



E 



ni 



k=i ^- pi+p2+...+pfc=n-i ^'I'^'a! •••Pfc' 



Proof: 

This will be done by recurrence. 
In fact, Aq = H{uq), 

then Ui = H(uQ)ds = H(uQ)t, and ui(x) = H(uo(x)) 
Jo 

and Uk{x,t) = Uk{x) 



Let's assume that \/k < n A^-i = uAx) 



{k-l)\ 



From Theorem 1, we have: = E IJ E -^''^^ (^o) (^pi 

fc=l ■ Pl+P2 + ---+Pfe=fl 



k\ 



Then: 



A 1 



fc=l Pi+P2+..-+Pfe=n 

E ry 

fc=l Pi+P2+--.+Pfc=n 

- E- E 



' ^P2 ' 



n\ 



fPl fP2 fPk 



PI' P2'-' P.y - p^lp^l-pfc! 



Pfc 



We have: 



n! 



n -t 

eA e - , , 

k=l ^' pi+p2+...+Pfe=n Pl'Ps! •••Pfc! 



^^''^ Kj iWpi '"p2 "k; ' 



and, M„+i = / Ands = Un+iix) / —ds 
Jo Jo nl 



this leads to the result. 



These two formulas show that Un closely depends on the derivatives H^''^ ("s"^) ' 
diiQ 

consequently on — — . On the other side the form of m„(x, t) = — can help us 

ox n'. 

oo 

to set some conjectures about the convergence of the serie Ukjx, t). 

k=0 



We have: 



Un+l 



Un+l{x) 



Unix 



t 



n + 1 



< 



Un+l{x) 



Unix 



\t\ 



Using the well known rule of d'Alembert, one can see that the convergence of 

°° ( d \ d 

uAx.t) depends on if'-'^^ I —Uq \ and —Uq . 
t^o V-9a; ) dx 



2 Characteristics method 



We still consider the problem in ([T]). 

By deriving formally in space (P^) and taking v 



9u , . 

— , we have the conservation 
ox 



4 



law equation 

7)^ -I- 

dx 



Vt + i-[H (v)] = (3) 



with the initial condition t>(x, 0) = Vo{x) = — 

ox 

Assuming that v is derivable in space, we then have: 

vt + H' (v) v, = (4) 
Using a(t>) = H'{v), we have: Vt + a (v) Vx = 

Let's consider in the plan {x,t) the characteristic strips x ^ x{t) solutions of the 
equation 



dx 

X(0) = Xq 

Writing w(t) = v {x(t), t), we have: 



— - a{v{x{t),t)) 



d . ^ d ,, , , / dx dv\ ,, , , 
-Mt) = -.((x,t),t) =(-. + ^^j((-,t),t) 

= {^t + a^ {{x,t),t)= 

This leads to: w(t) = wq] i.e. v is constant on the the characteristic strips. 
dx 

From ([5]) we have ~ (^o (^o)); then the characteristic strip from xq is a straight 
line with equation given by: 

x{t) = a {vq (xo)) t + xo (6) 




If a (uo (.)) is non decreasing, we obtain regular solutions , else one can observe losts 
of regularity (see figH]), i.e. the solution lost in regularity in time at T. 
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The time T* = inf{T} is called critical time of the equation ([3]) and is calculated by 



1 



. ( ^, duQ\ min (a(t;o)) (x) 

mm H o — — [x] x&R 



(7) 



dx 



du(] 



We can see that H' o — is one of the terms H^^"^ o — in the Adomian polynomials 

ox ox 

oo 

This leads us to say that the radius of convergence in time of the series ^ Uk{x,t) 
links closely to the critical time T*. 



3 Applications 



3.1 Example 1 

Let's consider the equation: 



1 2 

+ - (Ux) = 



u{x, 0) = uo{x) = — 

The associate conservation law problem is the well-known Burger's equation [TT] : 



Vt + 







v{x,0) = Vo{x) = —2x 



1 



We then have H{v) = -v^. Thus H'{v) = v = a(f), and: 
„, 1 1 



min (a(t;o))'(a;) 



min (vn) (x) 2 



Consider Theorem 1, the ADM gives: 



Mo 


= — 




= ^Kr = 


U\ 


= - J Aodt 








= -At'^x^ 


Un 


= -x2(2t)" 
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So: 

oo 

u{t,x) = -x^Y.(2tY- 

n=0 

For \t\ < 1/2 we have: 

u{t, x) 



x^ 



1 - It 



x^ 



One can easily verify that after t = 1/2, u{t,x) = — — remains a solution of our 

equation. But this is not the physically true. As we can see, u{t,x) isn't regular at 
t = 1/2. Recalling the fact that this equation model frontier evolution (for example 
in mathematical morphology [IZIIH]), we know that this irregularity spreads forward 
for if: > 1/2. At this stage, only generalized or entropic solutions survived after 

T* = -. We get these solutions by using Kruskov's formula (see for example f6]). 

Numerically these solutions can be computed with Godunov or Hamilton- Jacobi 
schemes. (see [T2| |9]) 



3.2 Example 2 

' ut + \/l + ul = 

, m(x, 0) = Uo{x) 
Case 1: Uo{x) = ax + b, a, 6 G IR. 

We can notice that the characteristic strips, in this case, are parallel lines. So this 
equation has a regular solution at all time. The viscosity solutions theory (see 
[TUl [T3] ) stipulates that the solutions are parallel lines to the initial condition. One 
can verify then H' (mq) is constant. Then T* = +oo. 
With the ADM, we have: 



Uo 



ax + b 



An 



Ui 



ty/T+a^ 



U2 



Ur 



Thus: 



Vn > 2 



u(t, x) = ax + b + tyl + a? 



This solution is exactly the viscosity solution given by the well-known Lax formula 



Case 2: Uq{x) = sinx 
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We then have H{v) = ^/l + v^. 
Thus H'(v) 



a[v\ 



and: T* 



min (a(fo))'(x) 



d 



mm 



^^0 



xglR dx \ Jl _|_ 



(x) 



d_ I Vq 
It is easy to see that: 



coss 



— smx 



dx \ Vl + cos^ X J (1 + cos^ x) y/1 + cos^ x 



— smx 



mm 



xGlR (1 + COS^ x) Vl + COS^ X 

Thus T* = 1. 
The ADM gives: 



— sm X 



[1 + cos^ x) Vl + cos^ X 



Ui 
M4 



smx 



+ COS^ X 



cos^ X sin a; 



2 1 + cos^ a; 

6(1 + COs2x)7/2 



(2 cos^ a; sin^ x + cos"^ a; sin^ x — cos'^ a; — 2 cos® x — cos^ a;) 



— ; — — (10 cos^ X sin^ X + 13 cos'^ x sin x — 4 cos^ x sin^ x 

24(1 + cos2x)5 ^ 

+24 cos® X sin x + 8 cos^ x sin x + 7 cos^° x sin x — 6 cos^ x sin^ x 

—6 cos^ X sin^ x + 3 cos^ x sin x) 





space: x 
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Figure 3: ^Uk{x,t) for (a) 0.5 < t < 1 and (b) l<t<l{) 

fe=0 



Figure 4: (a) Numerical solution with Hamilton- J acobi scheme and (h) Ukjx, t) at different 
times ( upwards : t = 0, 0.1, 0.25, 0.5, 0.75, 1, 1.25, 1.5, 1.75, 2, 2.25, 2.5;. 

One can see as shown on figures [3] and H] (a)-(b) that for < t < T* (with T* = 1 
the ADM gives an accurate approximation of the solution. 

This confirms that the ADM, for the case of Hamihon- Jacobi equations, gives a good 
approximation of the solution for only t <T* even with few terms Un- Nevertheless 
at t > T* the computation of several terms of the series doesn't ameliorate the 
solution. 



4 Conclusion 

We show in this paper that, for Hamilton- Jacobi equations and generally hyperbolic 
equations the ADM gives accurate solution or a good approximation of the solution 
only for t < T* . We deduce from there that ADM doesn't gives global solution in 
time for this type of equation. 
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